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In flooded and non-flooded impounded forests of Black mangrove {Avicennia germinans), 
the community structure of the ammonia-oxidizing betaproteobacteria (p-AOB) differed 
among distinct mangrove vegetation cover types and hydrological regimes. This had been 
explained by a differential response of lineages of p-AOB to the prevailing soil conditions 
that included increased levels of moisture and ammonium. To test this hypothesis, slurries 
of soils collected from a flooded and a non-flooded impoundment were subjected to 
enhanced levels of ammonium in the absence and presence of additional shaking. After 
a period of 6 days, the community composition of the p-AOB based on the 16S rRNA 
gene was determined and compared with the original community structures. Regardless 
of the incubation conditions and the origin of the samples, sequences belonging to the 
Nitrosomonas aestuarii lineage became increasingly dominant, whereas the number of 
sequences of the lineages of Nitrosospira (i.e., Cluster 1) and Nitrosomonas sp. Nm143 
declined. Changes in community structure were related to changes in community sizes 
determined by quantitative PCR based on the amoA gene. The amoA gene copy numbers 
of p-AOB were compared to those of the ammonia-oxidizing archaea (AOA). Gene copy 
numbers of the bacteria increased irrespective of incubation conditions, but the numbers 
of archaea declined in the continuously shaken cultures. This observation is discussed 
in relation to the distribution of the p-AOB lineages in the impounded Black mangrove 
forests. 
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INTRODUCTION 

Aerobic ammonia-oxidizing archaea and bacteria play an impor- 
tant role in the global nitrogen cycle by converting ammonia 
to nitrite (Kowalchuk and Stephen, 2001; Schleper and Nicol, 
2010). Nitrite can then be further oxidized to nitrate by nitrite- 
oxidizing bacteria or used as an electron acceptor in a large 
number of microbial and chemical redox reactions. Studies apply- 
ing the phylogenetic 16S rRNA gene and the amoA gene that 
codes for the ot-subunit of the enzyme ammonia monooxyge- 
nase in aerobic ammonia-oxidizing archaea and bacteria showed 
that the ammonia-oxidizing thaumarchaea and betaproteobac- 
teria (P-AOB) are widely distributed (Kowalchuk and Stephen, 
2001; Schleper and Nicol, 2010) and occur in many habitats, 
including soils of mangrove forests (e.g., Wickramasinghe et al., 
2009.) Mangroves are tree-dominated intertidal wetlands along 
tropical and subtropical coastlines with a specialized flora adapted 
to waterlogged and saline conditions. 

In a previous study (Laanbroek et al., 2012) we investigated 
the effects of vegetation cover type and summer flooding on the 
process of ammonia oxidation in soils collected from mangrove 



stands that were dominated by Black mangrove {Avicennia ger- 
minans) of differing statures that occurred in impoundments 
with dissimilar hydrologic regimes. The vegetation cover type did 
not significantly affect the potential ammonia-oxidizing activity 
(PAA) and a possible effect of summer flooding was offset by the 
effect of an extremely dry winter yielding high pore water salin- 
ities. In the same study (Laanbroek et al., 2012) we investigated 
also the distribution of different lineages of ammonia-oxidizing 
betaproteobacteria (P-AOB) in soils from different mangrove veg- 
etation cover types. Based on the 16S rRNA gene, the majority 
of the P-AOB in soils from these impoundments was related 
to the lineages of Nitrosomonas aestuarii/Nitrosomonas marina, 
Nitrosospira (Cluster 1) and Nitrosomonas Nml43, irrespective of 
the hydrologic regime. Sequences related to the N. aestuarii/N. 
marina lineage were mostly encountered at drier locations in 
the impoundments. In contrast, sequences associated with the 
Nitrosomonas Nml43 lineage were mostly observed at the wet- 
ter sites. Since higher moisture levels are generally associated 
with more reduced conditions in the soils, it can be argued that 
members of the Nitrosomonas Nml43 lineage are better adapted 
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to such reduced conditions. Sequences related to Nitrosospira 
Cluster 1 were mostly found at places with low nutrient availabil- 
ity as had been observed before in a freshwater marsh (Laanbroek 
and Speksnijder, 2008). Hence, the lineages of P-AOB were differ- 
ently affected by a number of environmental factors. 

To test the hypothesis that ammonium availability and redox 
conditions impacted the characteristics of the f5-AOB com- 
munity in the soils dominated by Black mangrove, changes 
in the composition of the communities were determined in 
ammonium-enriched, shaken, or non-shaken soil suspensions 
that were incubated for 6 days at room temperature. It was 
assumed that shaking resulted in more oxic conditions for the 
aerobic ammonia-oxidizing cells. Coci et al. (2005) had previ- 
ously shown that the species composition of P-AOB communities 
can respond quickly to changes in environmental conditions 
and that oxic conditions are a prerequisite for changes in the 
community structure. Based on their distribution in the Black 
mangrove soils, we expected that gene numbers of Nitrosomonas 
lineages would increase in the nutrient-enriched soil suspen- 
sions at the cost of Nitrosospira Cluster 1, irrespective of shaken 
or non-shaken conditions, and that shaking of the suspensions 
will give the N. aestuarii/N. marina lineage extra advantage over 
the Nitrosomonas Nml43 lineage. Changes in community struc- 
ture of the P-AOB due to 6 days prolonged incubation were 
related to changes in community sizes determined by quantita- 
tive PCR based on the amoA gene. The amoA gene copy numbers 
of P-AOB were compared to those of the ammonia-oxidizing 
archaea (AOA). 

MATERIALS AND METHODS 
STUDY LOCATIONS 

The two mangrove-dominated impoundments that we studied, 
and described in detail in Laanbroek et al. (2012), are located in 
the Indian River Lagoon, a coastal sub-estuary, located on North 
Hutchinson Island, St. Lucie County, Florida, USA. For the exper- 
iments and analyses described in the present paper, samples were 
collected in March 2009 and March 2010 just after a new man- 
agement regime of spring and summer flooding was initiated for 
purposes of controlling the numbers of salt-marsh mosquitoes 
and biting midges in one of the impoundments (Impoundment 
24). Hence, soils collected in March 2009 and March 2010 expe- 
rienced different hydrologic regimes in Impoundment 24 during 
the preceding spring and summer period resulting in increased 
days of total soil submergence (Verhoeven et al., under review). 
For comparison, soil samples were also collected from the adja- 
cent Impoundment 23 that was not subjected to the new man- 
agement. Soil samples used in this study were collected from 
three different vegetation cover types, each dominated by Black 
mangrove (Avicennia germinans). The three cover types (dwarf, 
sparse, and dense) are described in Laanbroek et al. (2012) and 
Verhoeven et al. (under review). Within each impoundment, soils 
from five locations in each of the three cover types were collected 
resulting in a total of 30 locations sampled. 

SOIL COLLECTIONS 

In March 2009 and March 2010, soil cores (3.9 cm diameter and 
10 cm long) were collected at each of the 30 sites. Samples were 



collected with an aluminum tube that was sharpened at one end. 
The cores were immediately sealed at both ends with rubber stop- 
pers and transported to the laboratory. One of the cores was used 
for removal of pore water that was analyzed for salinity and pH. 
The top 5-cm of two other cores were combined and thoroughly 
mixed by hand and subsequently sub-sampled for the determi- 
nation of bulk density and moisture content, for measurements 
of potential rates of nitrification and for the analysis of the com- 
munity structure of P-AOB. Samples for the genetic analyses were 
freeze-dried immediately after mixing of the sampling cores and 
stored until further processing. 

SLURRY EXPERIMENTS 

PAAs were determined in slurries of 20 g fresh weight soil mixed 
with 50 ml of mineral medium containing ammonium at a final 
concentration of 1 mM (see Laanbroek et al., 2012 for details). 
Central to this PAA method is the determination of the accumula- 
tion rate of nitrite plus nitrate being together the sum of ammonia 
oxidized by the microorganisms present in the soil samples. After 
the initial measurements of 9h under optimal conditions (i.e., 
shaken) at room temperature, the slurries were left at room tem- 
perature for 6 days, either non-shaken (2009 samples) or shaken 
(2010 samples). At day 6, the PAA measurements were repeated 
after ImM ammonium had again been added, while all slurries 
were aerated again; hence, both in 2009 and 2010. After the last 
sampling on day 6, soil particles in the slurries that had settled out 
for 1 h were removed and subsequently freeze-dried for genetic 
analyses. In total, 30 slurry samples were measured each year, 
i.e., slurries from quintuplet soil samples from three different 
vegetation cover types in two different impoundments. 

DNA ISOLATION 

DNA isolation, PCR, and construction of clone libraries were sim- 
ilar to the procedures described for the freshly collected soil sam- 
ples as described before (Laanbroek et al, 2012). From each of the 
slurry samples that were freeze-dried after the last measurement 
on day 6, 0.5 g was homogenized in 1 ml cetyltrimethylammo- 
nium bromide (CTAB) buffer (Zhou et al, 1996) in the presence 
of 0.5 g sterilized zirconia-silica beads (diameter 0.1mm), sub- 
jected to disruption by bead-beating at a 5.0-m/s rotation for 60 s 
and subsequently incubated for 30min at 37° C in the presence 
of 5 |xl proteinase K (20 mg/ml). After vortexing for 15 min, the 
samples were supplemented with 150 (xl of a 20% SDS solution 
and incubated for 1 h at 65°C in a Thermoblock and again vor- 
texed every 15-20 min. After centrifugation at 10,000 x g for 
10 min, 600 jxl of supernatant was collected in 2-ml screw-cap 
tubes. The rest of the sample was re-extracted with 450 [i\ CTAB 
buffer and 50 jxl of a 20% SDS solution, vortexed for 10 s, incu- 
bated for 10 min at 65°C, and centrifuged at 6,000 x g for 10 min. 
Again, 600 uJ was collected, added to the previously extracted 
supernatant, mixed with 1 ml phenol-chloroform-isoamyl alco- 
hol solution (25:24:1, vol/vol/vol), and centrifuged at 6,000 x g 
for 10 min. One milliliter of supernatant was collected and placed 
into a new screw-cup tube containing 700 (jlI isopropanol, and the 
tube was incubated for 1 h at 24° C. After 20 min of centrifugation 
at 15,000 x g, the isopropanol was decanted and the pellet was 
resuspended and washed with 1 ml 70% cold ethanol. This was 
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followed by 5min of centrifugation at 15,000 x g, decantation 
of the ethanol, drying of the pellet under vacuum centrifugation, 
and finally, resuspension in 100 [il water (Sigma). 

PCR AND CONSTRUCTION OF CLONE LIBRARIES BASED ON THE 16s 
rRNAGENE 

Extracted DNA was amplified by nested procedure using two 16s 
rRNA gene primer sets specific for the majority of the betapro- 
teobacterial AOB, i.e., the pAM0161f and f5AMO1301r primer 
set of McCaig et al. (1994) and the CTC)189f and CTO654r 
primer set of Kowalchuk et al. (1997). One hundred ng of purified 
DNA was used as template for a 50-u,l PCR mixture contain- 
ing 1 x Mg-free buffer (Invitrogen Corp., Carlsbad, CA), 0.5 u,M 
of each primer, 200 |iM of each deoxynucleotide triphosphate, 
1.75 mM MgCi2, 400 ng/(xl bovine serum albumin, 1.25 U GoTaq 
Hot Start Polymerase (Promega). The thermocycling program for 
both steps consisted of 2 min of denaturation at 95°C followed by 
35 cycles of 30 s of denaturation at 95° C; 30 s of specific annealing 
at 59°C (PAMO primer set) or at 57°C (CTO primer set) and 45 s 
of elongation at 72°C; 5 min of final elongation was performed for 
all reactions. Nested amplifications of 25 cycles were performed 
with the primer set CTO189f and CTO654r on 1:100 dilutions of 
PCR products from the f3AMO primer set. All reactions were ver- 
ified by UV illumination of 1% agarose gels stained in a gel red 
or ethidium bromide solution. Polymerase chain reaction frag- 
ments were ligated into the pGEM T-vector system (Promega 
Corporation, Madison, WI, USA) and transformed into JM109 
competent E. coli cells (Promega) according to the manufactur- 
ers' instructions. Transformed colonies were screened for inserts 
of the correct size by PCR. Sanger sequencing using primer T7 
was performed by Macrogen, Amsterdam, Netherlands. 

SEQUENCE ANALYSES 

Sequences were aligned and checked for chimeras using 
Sequencher 4.1 (Gene Codes Corporation, Ann Arbor, MI). 
The chromatogrammes of the remaining sequences were visu- 
ally inspected and wrongly assigned bases were manually 
replaced. Low quality sequences were removed from the dataset. 
After removing sequences of insufficient length or quality, 
the remaining 1022 sequences were further analyzed. The 
aligned sequences were clustered in Operational Taxonomic 
Units (OTUs) with the Cluster program of MOTHUR software 
version 1.18.1 (Schloss et al., 2009). In addition, MOTHUR 
was also used for establishing rarefaction curves and group 
(i.e., OTU-) representatives. OTUs at the level of 97% mutual 
similarity in their 16S rRNA gene sequences (418 bp) were 
chosen for further analyses. Sequence identification of the OTU- 
representatives was done by the BLASTN facility from the 
National Center for Biotechnology Information. (http://www. 
ncbi.nlm.nih.gov/). The partial 16SrRNA gene sequences were 
submitted to the GenBank database under accession numbers 
JQ725556— JQ726359 (fresh soils) and KF7 1 9449— KF7 19951 
(incubated slurries). 

Analysis of similarity (ANOSIM) between communities of f5- 
AOB obtained from different impoundments, years, mangrove 
cover types and incubations times of the slurries was done 
with the PRIMER software version 5.2.9 (Primer-E, Plymouth, 



UK). The ANOSIM analysis was based on Bray-Curtis distance 
matrices. 

QUANTITATIVE PCR OF ARCHAEAL AND BACTERIAL amoA GENES 

Gene copy numbers of archaeal and bacterial amoA genes were 
determined in separate reactions with the primer sets ArchamoA- 
1F and ArchamoA-2R (Francis et al., 2005) and amoA-lF and 
amoA-2R (Rotthauwe et al., 1997), respectively. The reaction 
mixture consisted of 10 \iL Rotor-Gene SYBR Green PCR kit 
(Qiagen), 1 (xL of each primer (5 pmol (JtL -1 ), 10u,g bovine 
serum albumin (100 mg ml , New England Biolabs Inc.) and 
5 |xL template DNA (approximately 50 ng total soil DNA); final 
volume was 20 ul. The PCR conditions for the archaeal amoA 
gene were 15 min at 95°C, followed by 40 cycles of 45 s at 95°C, 
45 s at 55°C and 45 s at 72°C, and finally 5 min at 72°C. The PCR 
conditions for the bacterial amoA gene were 5 min at 95°C, fol- 
lowed by 45 cycles of 45 s at 95° C, 45 s at 58°C and 45 s at 72°C, 
and finally 5 min at 72° C. Standard curves for quantification of 
bacterial and archaeal amoA genes were generated fromlO-fold 
serial dilutions (10 2 -10 8 copies |iL _1 ) of plasmid DNA contain- 
ing either a bacterial amoA fragment of 635bp that was obtained 
from a pure culture of Nitrosomonas europaea ATCC 19718 or 
an archaeal amoA fragment from clone 29C_47 [accession num- 
ber JQ404089, (Daebeler et al, 2012)]. The detection limits for 
the qPCR assay were 2.83 x 10 3 and 9.81 x 10 2 gene copies/g of 
dry soil for the archaeal and bacterial amoA genes, respectively, 
which corresponds to three copies per reaction. qPCR was per- 
formed in a Corbett Research Rotorgene 3000 Cycler, (Qiagen, 
Netherlands). 

STATISTICAL ANALYSES 

In contrast to the log-transformed bacterial amoA gene copy 
numbers, the log-transformed archaeal amoA gene copy num- 
bers were not normally distributed. Therefore, nonparametric 
Spearman rank order correlation and Kruskal-Wallis analysis of 
variation were applied for both archaeal and bacterial amoA gene 
copy numbers. The statistical analyses were performed with the 
STATISTICA software version 12 (StatSoft Inc., Tulsa, OK). 

RESULTS 

POTENTIAL AMMONIA-OXIDIZING ACTIVITIES 

PAAs in the soil slurries increased significantly (p = 0.0039 and 
p < 0.0001) in 2009 and 2010 with prolonged incubation for 6 
days in the laboratory (Figure 1). The increase in activity was 
most pronounced in 2010 when shaking of the ammonium- 
enriched slurries continued during the whole incubation period. 
In 2009, shaking had been discontinued between the potential 
activity measurements on day 0 and 6. 

OTUs OF p-AOB PRESENT IN MANGROVE SOIL SAMPLES 

Based on 97% mutual similarity between the 16S rRNA genes 
(418 bp), the 1067 qualified sequences could be classified into 35 
different OTUs. Of these 35 OTUs, 17 contained only a single- 
ton, and 10 others comprised only 2-4 sequences. The remaining 
eight OTU's contained at least eight sequences and were used 
for the analyses of community composition. These eight OTUs 
comprised 320 and 220 sequences obtained from the original 



www.frontiersin.org 



November 2013 | Volume 4 | Article 343 | 3 



Laanbroek et al. 



Response of ammonia-oxidizing bacteria to environmental changes 



mangrove soils in 2009 and 2010, respectively, and 221 and 261 
sequences from the incubated slurries at the end of the incuba- 
tion experiments in 2009 and 2010, respectively. A BLAST analysis 
of the eight OTU's representatives yielded a limited number of 



lineages of P-AOB (Table 1). Except for OTU05, which sequences 
were 100% similar to a sequence obtained from a nitrifying 
bioreactor, all other OTUs contained sequences that were mostly 
related to sequences originating from estuarine, marine or saline 
environments. 

INCUBATION-INDUCED CHANGES IN COMMUNITY COMPOSITION IN 
SLURRIES 

Incubation of the soil slurries in the absence of shaking led to a 
significant change in community composition (ANOSIM: R = 
0.084, p = 0.050). The average dissimilarity between the com- 
munities at the start and the end of the incubation period was 
59%. OTU01 contributed most (i.e., 39%) to this dissimilar- 
ity followed by OTU02 and OTU04 with 24 and 15%, respec- 
tively (Figure 2, supplementary Table 1). The relative presence of 
sequences belonging to OTU01 (related to N. aestuarii) increased 
from 46 to 69%. With an increase of 6-22% of the total sequences, 
the increase in OTU01 was most pronounced in the samples from 
locations with the sparse vegetation cover types. In the same 
period of 6 days, the relative occurrence of sequences belonging 
to OTU02 (related to Nitrosospira Cluster 1) decreased from 22 
to 9% of the total number of sequences. This relative decrease 
in presence was observed for all vegetation cover types. The rel- 
ative changes in occurrence of sequences belonging to the other 
OTUs were small except that sequences belonging to OTU05 
(comprising the N. europaea lineage) and OTU07 (consisting of 
sequences mostly related to the N. aestuarii/N. marina lineage) 
that disappeared completely from the clone libraries. 

Incubation of the soil slurries in the presence of shak- 
ing led also to a significant change in community composi- 
tion (ANOSIM: R = 0.149, p = 0.014).The average dissimilarity 
between the communities at the start and the end of the incuba- 
tion period was 64%. OTU01 contributed again most (i.e., 35%) 
to this dissimilarity followed this time by OTU03 and OTU02 
with 24 and 17%, respectively (Figure 3, Supplementary Table 2). 
The relative presence of sequences belonging to OTU01 more 
than doubled from 26 to 54% of the total sequences. This was 
true for all three vegetation cover types, except for the dwarf veg- 
etation type where OTU01 comprised already almost half of the 
sequences before the incubation started. Sequences belonging to 
OTU02 almost disappeared from the clone libraries with only 
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FIGURE 1 | Average potential ammonia-oxidizing activities in 
mangrove soil slurries before and after prolonged incubation of 6 days 
with added ammonium in the absence (A) and presence (B) of 
continuously shaking of the slurries. Soils had been collected in 2009 
and 2010 in two different impoundments (23 and 24) from three different 
habitats characterized by stunted (blue), sparse (red), and dense (green) 
Avicennia germinans. Both in 2009 and in 2010, incubation increased the 
potential activities significantly (p = 0.0039 and p = 0.0000, respectively; 
non-parametric Kruskal-Wallis test with N = 60 and H = 8.344 and 
H = 27.876, respectively). No significant effects between years, 
impoundments, and vegetation cover types were observed. 



Table 1 | BLAST analysis of OTU's representatives encountered in Black mangrove soils collected from the Indian River lagoon at North 
Hutchinson Island near Fort Pierce, Florida. 



OTU# 




BLAST analysis of OTU's representatives 




Closest type strain 


% similarity 


Origin closest relative 


References 


01 


Nitrosomonas aestuarii 


97 


High altitude saline wetland 


Doradoretal., 2008 


02 


Nitrosospira tenuis 


96 


Deep water sponge 


Meyer and Kuever, 2008 


03 


Nitrosomonas sp. Nm143 


97 


Estuarine sediment 


Jin et al., 2011 


04 


Nitrosomonas aestuarii 


98 


Prawn farm sediment 


Ma and Wang, 2007 


05 


Nitrosomonas europaea 


100 


Nitrifying bioreactor 


Bae and Song, 2011 


06 


Nitrosomonas aestuarii 


97 


Estuarine sediment 


Magalhaes et al., 2007 


07 


Nitrosomonas aestuarii 


98 


Estuarine sediment 


Freitag et al., 2006 


08 


Nitrosomonas aestuarii 


96 


Coastal marine sediment 


Urakawa et al., 2006 



Only OTUs with numbers of sequences larger than 1 % of the total number of sequences. 
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FIGURE 2 | Effect of 6 days of ammonium enrichment on the sequence 
distribution of ammonia-oxidizing betaproteobacteria in non-shaken 
soil slurries. Sequences originated from different vegetation cover types of 
Avicennia germinans forests sampled in 2009. Operational Taxonomic Units 
(OTUs) are based on the partial 16S rRNA gene (440 bp). 



FIGURE 3 | Effect of 6 days of ammonium enrichment on the sequence 
distribution of ammonia-oxidizing betaproteobacteria in shaken soil 
slurries. Sequences originated from different vegetation cover types of 
Avicennia germinans forests sampled in 2010. Operational Taxonomic Units 
(OTUs) are based on the partial 16S rRNA gene (440 bp). 



a few sequences left in the dense vegetation cover type. With a 
decrease from 23 to 8% of the total, the relative occurrence of 
sequences belonging to OTU03 declined especially in the soil slur- 
ries originating from the dense vegetation cover type. The relative 
changes in occurrence of sequences belonging to the other OTUs 
were small. However, the few sequences belonging to OTU02 and 
OTU08 disappeared completely from the slurries. 

amoA GENE COPY NUMBERS BEFORE AND AFTER PROLONGED 
INCUBATION 

Averaged for the sampling groups, numbers of archaeal amoA 
genes were always larger than the numbers of bacterial amoA 
genes, irrespective of the year of sampling, the impoundment 
from which the samples were collected and the mangrove veg- 
etation cover type (Figure 4). For individual samples, the ratios 
between archaeal and bacterial gene copy numbers varied largely 
with the lowest (i.e., 0.9) and highest (6.5 x 10 4 ) ratios found in 
2009 in the sparse vegetation of impoundment 23 and impound- 
ment 24, respectively. Gene copy numbers of both AOA and 
P-AOB were significantly larger in 2010 compared to 2009 
(Table 2). The impoundment from which the soil samples had 
been collected had no significant effect on the gene copy numbers, 
although the ratio of archaeal to bacterial gene copy numbers 
were significantly higher in impoundment 24. On average, both 
archaeal and bacterial gene copy numbers were larger in the sparse 



mangrove vegetation compared to the other cover types; how- 
ever, there was no significant difference in the bacterial gene copy 
numbers between the sparse and the dense mangrove vegetation 
types. 

The archaeal and bacterial gene copy numbers were sig- 
nificantly correlated with each other (Supplementary Table 3: 
Spearman rank order correlation). Both correlated also signifi- 
cantly with extractable soil phosphate (positively) and with pore 
water pH and nitrate (negatively). The archaeal amoA gene copy 
numbers, but not the numbers of P-AOB was significantly and 
negatively correlated with pore water salinity. Finally, neither 
archaeal nor bacterial amoA gene copy numbers correlated with 
the PAA. 

Prolonged incubation for 6 days in the presence of added 
ammonium increased both the archaeal and the bacterial amoA 
gene copy numbers (Figure 4; 2009 data series). These increases 
were significant (Table 3). However, the copy numbers of the 
P-AOB increased more leading to a significantly lower AOA to 
P-AOB ratio. Prolonged incubation for 6 days in the presence 
of both added ammonium and shaking had a significant neg- 
ative effect on the archaeal amoA gene copy numbers, whereas 
the numbers of the bacterial amoA genes increased (Figure 4; 
2010 data series). Again, the ratio between archaeal and bacte- 
rial gene copy numbers declined significantly during the period 
of prolonged incubation. 
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AOA (3-AOB 



FIGURE 4 | Archaeal (left panels) and bacterial (right panels) amoA 
gene copy numbers (numbers per g dry soil) before (blue bars) 
and after (red bars) prolonged incubation of 6 days with added 
ammonium in the absence (2009 series) and presence (2010 series) 



of continuously shaking of the slurries. Soils had been collected in 
two different impoundments (23 and 24) from three different habitats 
characterized by stunted, sparse, and dense Avicennia germinans. For 
statistical data see Table 2. 



Table 2 | Effect of year of sampling, impoundment and type of Black 
mangrove vegetation on the archaeal and bacterial amoA gene copy 
numbers (AOA and AOB, respectively) and their mutual ratio 
according to a nonparametric Kruskal-Wallis test. 







N 


Degrees of 


P 


Remarks 








freedom 






Year 


log AOA 


55 


1 


0.0014 


2010 > 2009 




log AOB 


55 


1 


0.0959 


2010 > 2009 




log AOA/AOB 


52 


1 


0.6905 




Impoundment 


log AOA 


55 


1 


0.7738 






log AOB 


55 


1 


0.1340 






log AOA/AOB 


52 


1 


0.0197 


24 > 23 


Vegetation type 


log AOA 


55 


2 


0.0067 


SP > ST= DE 1 




log AOB 


55 


2 


0.0024 


SP = DE > ST 




log AOA/AOB 


52 


2 


0.0020 


SP = DE > ST 



1 ST stunted; S/? sparse; DE, dense. 



DISCUSSION 

The hypothesis that enrichment with ammonium and more oxic 
conditions would favor members of the N. aestuarii/N. marina 
lineage more than the other lineages present in the soil sam- 
ples was confirmed by the slurry experiments. During the slurry 
experiments, OTU01 related to theN. aestuarii/N. marinalineage, 
which was dominant in the original soil samples (Figures 1,2) 
became even more dominant at the cost of OTU02 belonging 
to Nitrosospira Cluster 1 and to a lesser extent at the expense of 
OTU03, which is related to the Nitrosomonas sp. Nml43 lineage. 



Table 3 | Effect of incubation on the numbers of ammonia-oxidizing 
archaea and bacteria and their mutual ratio according to a 
nonparametric Kruskal-Wallis test. 





Year of 


N 


Degrees of 


H 


P 


Remarks 




sampling 




freedom 








log AOA 


2009 


55 


1 


7.044 


0.0080 


day 6 > day 0 




2010 


59 


1 


15.603 


0.0001 


day 6 < day 0 


log AOB 


2009 


53 


1 


9.284 


0.0023 


day 6 > day 0 




2010 


59 


1 


0.568 


0.4439 


day 6 > day 0 


log AOA/AOB 


2009 


50 


1 


2.992 


0.0837 


day 6 < day 0 




2010 


58 


1 


5.126 


0.0236 


day 6 < day 0 



The changes observed in ammonium-enriched, non-shaken soil 
slurries were amplified strongly in ammonium-enriched, shaken 
incubations. The more oxic conditions of the soil suspensions 
due to shaking strengthened the effect of incubation on the rel- 
ative presence of the different lineages compared to the treatment 
without shaking. Apparently, the presence of both ammonium 
and more oxic conditions offered members of the N. aestuarii/N. 
marina lineage conditions for growth that enabled them to over- 
grow their competitors. The large decline with 15% in the relative 
occurrence of sequences belonging to OTU03 when exposed to 
shaking might reflect the adaptation of the Nitrosomonas sp. 
Nml43 to more reduced conditions as these likely prevail in 
the moist dense vegetation cover type in the impoundments 
(Verhoeven et al, under review). Such a large decline in the 
relative occurrence of sequences belonging to OTU03 was not 
observed in the non-shaken soil slurries. During incubation of 
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the soil slurries, the relative decline in the number of sequences 
related to the Nitrosomonas sp. Nml43 lineage was smaller than 
the decline in the number of sequences belonging to Nitrosospira 
Cluster 1. This could be due to a yet unknown characteristic that 
enables Nitrosomonas sp. Nml43 to counterbalance the dominant 
N. aestuarii/N. marina lineage better or for a longer period than 
the Nitrosospira lineage. 

Incubation of the slurries resulted in an increase in the PAAand 
in the community composition of f5-AOB over a period of 6 
days. The increased activity was apparently related to the observed 
increases in bacterial amoA gene copy numbers in both years. 
Similar to the bacterial gene numbers, the archaeal amoA gene 
copies also increased in numbers during the prolonged 6 days 
incubation period, but only when ammonium had been added 
in the absence of prolonged shaking. When the soil slurries were 
continuously shaken during the incubation period, archaeal gene 
copy numbers declined. Hence, AOA may be less resistant to more 
oxic conditions as caused by continuously shaking of the soil sus- 
pensions. Data on ratios between AOA and P-AOB gene copy 
numbers in mangrove forest soils are relative scarce and noth- 
ing is known about the effect of redox conditions on this ratio. 
The numbers of bacterial amoA genes in a Kandelia obovata- 
dominated mangrove forest in the Mai Po Nature Reserve in the 
New Territories of Hong Kong outnumbered those of archaeal 
amoA genes by a factor of 5-12 (Li et al, 2011). In mangrove 
soil microcosms from the same area, addition of addition of 
ammonium led to an increase in amoA gene copy numbers of 
both AOA and f5-AOB, but the latter remained numerically dom- 
inant (Li and Gu, 2013). Studies with respect to the presence 
of ammonia-oxidizing gammaproteobacteria in mangrove for- 
est soils are limited. Sequences belonging to this bacterial group 
could not be demonstrated in this habitat ( Wickramasinghe et al., 
2009). 

Aerobic degradation of labile materials near the surface of 
mangrove sediments happens usually so fast that oxygen rarely 
penetrates more than 2 mm into the sediment (Kristensen et al., 
1994). Redox potentials of Avicennia-dominated sediments at the 
coast of Queensland, Australia, varied from 0 till —200 mV at 
the surface and reached values from —100 till — 300 mV at 5 cm 
depths (Clark et al., 1998). The Avicennia-dominated forest soils 
that we use in our incubation experiments had been exposed to 
atmospheric oxygen for several months prior to sampling. This 
was particularly true for the dwarf and sparse vegetation cover 
types. Nevertheless, it is likely that a part of the aerobic ammonia- 
oxidizing community in the upper 5 cm of the soil that had been 
collected was inactive at the time of sampling due to oxygen 
limitation. 

What do the observations tell us about the distribution of 
lineages of fi-AOB in the mangrove soils from which the sam- 
ples were collected? Sequences belonging to the N. aestuarii/N. 
marina lineage were most numerous among the f5-AOB encoun- 
tered in the Black mangrove soils from the impoundments along 
the Indian River lagoon. They accounted for nearly 50 percent of 
the P-AOB community. Hence, it maybe concluded that the con- 
ditions in the impoundments are especially suitable for the N. aes- 
tuarii/N. marina lineage. Next to the conditions indicated above, 



salinity may be an additional factor favoring the N. aestuarii/N. 
marina lineage. The share of sequences belonging to OTU01 was 
higher in 2009 compared to 2010 (Tables 2, 3), which might have 
been due to average higher salinities measured in 2009 compared 
to 2010. However, the slurry experiments were performed at their 
own field salinity, which means in ranges of 52-88% and 28-54% 
in 2009 and in 2010 respectively, which likely exclude salinity as 
factor that induced changes in the communities of f$-AOB during 
slurry incubations. 

The Black mangrove forests in the impoundments along 
the Indian River lagoon are rather atypical for this type of 
ecosystem in the sense that tide is largely reduced. In addition, 
Impoundment 24 had been submerged during the summer sea- 
son of 2009, which resulted in higher ammonium concentrations 
in the pore water (Verhoeven et al., under review). Especially 
the dwarf and sparse black mangrove vegetation cover types 
have been exposed directly to the atmosphere for several months 
during the autumn and winter period before the samples were 
collected. Under these exposed conditions with increased ammo- 
nium concentrations, members of the N. aestuarii/N. marina 
lineage become apparently the dominant f5-AOB. 

The lineages of P-AOB as detected in the Black mangrove 
soils are commonly observed in estuaries (Francis et al., 2003; 
Bernhard et al, 2005, 2007; Freitag et al, 2006; Ward et al, 
2007; Mosier and Francis, 2008; Sahan and Muyzer, 2008; Jin 
et al., 2011; Wankel et al., 2011). Their mutual distribution is 
strongly influenced by salinity, with the Nitrosomonas lineages 
more at the land side and the Nitrosospira lineages more at the sea 
side of the estuaries. However, as already concluded by Francis 
et al. (2003), no single physical or chemical parameter entirely 
explains the pattern of diversity along the estuary, suggesting 
that a complex combination of environmental factors (e.g., oxy- 
gen, temperature, nitrate, and ammonium concentrations) may 
shape the overall level of AOB diversity in a dynamic environ- 
ment. In our paper, we have indeed shown that factors such 
as ammonium availability and redox conditions play indeed a 
role in structuring communities of P-AOB as well. The ecologi- 
cal consequences of a change in community composition of the 
P-AOB require further study. A change in community composi- 
tion of P-AOB in the Scheldt estuary due to increasing salinity 
(De Bie et al., 2001) had in fact no effect on the rate of ammo- 
nia oxidation as this was only determined by the availability 
of ammonium and oxygen, which were both limiting microbial 
activity. However, a replacement of lineages of f5-AOB was essen- 
tial for the maintenance of the ecological function of ammonia 
oxidation at changing environmental conditions. The same phe- 
nomenon may apply for the studied impounded Black mangrove 
soils. 
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